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The impact of meridional circulation on stellar butterfly diagrams 
and polar caps 
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ABSTRACT 

Observations of rapidly rotating solar-like stars show a significant mixture of opposite-polarity 
magnetic fields within their polar regions. To explain these observations, models describing 
the surface transport of magnetic flux demand the presence of fast meridional flows. Here, 
we link sub-surface and surface magnetic flux transport simulations to investigate (i) the 
impact of meridional circulations with peak velocities of < 125 m- s _1 on the latitudinal 
eruption pattern of magnetic flux tubes and (ii) the influence of the resulting butterfly dia- 
grams on polar magnetic field properties. Prior to their eruption, magnetic flux tubes with low 
field strengths and initial cross sections below ~ 300 km experience an enhanced poleward 
deflection through meridional flows (assumed to be poleward at the top of the convection 
zone and equatorward at the bottom). In particular flux tubes which originate between low 
and intermediate latitudes within the convective overshoot region are strongly affected. This 
latitude-dependent poleward deflection of erupting magnetic flux renders the wings of stel- 
lar butterfly diagrams distinctively convex. The subsequent evolution of the surface magnetic 
field shows that the increased number of newly emerging bipoles at higher latitudes promotes 
the intermingling of opposite polarities of polar magnetic fields. The associated magnetic flux 
densities are about 20% higher than in the case disregarding the pre-emptive deflection, which 
eases the necessity for fast meridional flows predicted by previous investigations. In order to 
reproduce the observed polar field properties, the rate of the meridional circulation has to be 
on the order of 100 m • s _1 , and the latitudinal range from which magnetic flux tubes originate 
at the base of the convective zone ( < 50°) must be larger than in the solar case ( < 35°). 

Key words: stars: magnetic fields - stars: activity - stars: interior - stars: rotation - stars: 
spots - stars: imaging 



1 INTRODUCTION 

On the Sun, dark spots are exclusively found within an equato- 
rial belt between about ±40° latitude. In contrast, stellar surface 
bright ness maps, secur ed with the technique of Doppler imaging 
ICoUier Cameroal200lL and references therein), show that other 
cool stars frequently have large high-latitude and pol ar spots, of- 
ten in conjunction with low-latitude features as well l Strassmeier 
2002, and references therein). Due to specific requirements of the 
observing technique the targets so far are stars rotating more rapidly 
than the Sun. Although the time base of surface maps acquired for 
individual stars is yet not sufficiently long to conclusively discern 
long-term activity properties, their behaviour seems to digress from 
the solar 1 1-year spot/22-year magnetic cycles and polar field prop- 
erties. In fact, many of the younger and more active stars show no 
apparent cycle insofar as chromospheric indicators such as Ca II H 
& K can be used as a proxy for magnetic activity I Baliunas et al. 
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1995; Donahue 1996). There is, as yet, also no example of a star 
undergoing a 'Maunder minimum'. Different activity and cycle 
signatures may however be present or more pronounced, such as 
apparent preferred longitudes in the spot distribution or the 'flip- 
flop' phenom enon l le tsu et alll994t lBerdvugina 2004; Moss 2004; 
iKorhonen & ElstneJ2005l) ! 

An important characteristic difference to the solar-like 
magnetic field distribution is the significant mixture of mag- 
netic flux of opposite polarity within polar regions, as ob- 
served with the Zeeman-Doppler Imaging technique (Semel 1989; 
Donati & Brown 1997). In particular the strong flux intermin- 
gling in the polar regions is in contrast to the Sun JDonati et all 
2003), where the high-latitude magnetic field is essentially unipo- 
lar throughout the majority of the activity cycle. On rapidly rotating 
stars, the associated polar magnetic flux densities are sufficiently 
high and persistent to cause dark pol ar caps lasting ove r a large 
number of stellar rotation periods (e.g. lleffers et al . 2005). 

The magnetic activity signatures of cool stars are ascribed 
to the emergence of magnetic flux generated by sub-surface dy- 
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namo mechanisms. Scenarios for dynamo operation are based on 
the magneto-hydrodynamic interaction between convective mo- 
tions and (differential) rotation. Yet there is no complete theory, 
which unifies the amplification, storage, transport, eruption, and 
(possibly cyclic) re-generation of magnetic fl ux consistently; for a n 
extensive review on stellar dynamo theory see Ossendrijver 1 2003). 
Dynamo mechanisms inside the convective envelope leave charac- 
teristic imprints on the observable activity signatures in the stel- 
lar atmosphere, which provide constraints for the underlying pro- 
cesses. 

In the following, we focus on the transport of magnetic 
flux both bel ow and on the stellar surface. Magnetic flux tube 
models (e.g. Moreno-Insertis 1986; Choudhuri & Gilman 1987; 
lD'Silva& Choudhuri 19931: iFan et alJll994ISchussler et alll994 



1996) describe the evolution of magnetic flux, concentrated in 
strands of magnetic field lines, inside the convection zone until 
their eruption on the stellar surface. Although they exclude as- 
pects concerning the (cyclic) generation of magnetic fields, they 
successfully reproduce characteristic properties of emerging bi- 
polar spot groups like, in the case of the Sun, their latitude of 
emergence, their relative velocities and topologies, the asymme- 
tries bet ween the prece ding- and following spot group, and Joy's 
law (e.g. iD'Silva & Choudhurilll993l: ICaligari et al J Il995l> . Fur- 
ther applications of the flux eruption model comprise cool stars 
with different rotation rates, stellar masses, and evolutionary stages 
as well as compone nts of close binary systems l Granz er et alJ 
l2000UHolzwarth & SchussleJl200ll E)03: Holzwarth 2004J). Fre- 
quently observed starspots at higher latitudes of rapid rotators, 
for example, can be explained by the poleward deflection of ris- 
ing fl ux tubes prior to the ir eruption on the stellar surface (e.g. 
Schiissler & Solanki 1992). A persistent magnetic flux eruption 
at high latitudes implies however the availability (and possibly 
generation) of large amounts of magnetic flux over a similar lat- 
itudinal range inside the stellar convection zone, which would 
be in dissent with the solar case, where the production of mag- 
netic flux is anticipated to be most efficient at low latitudes (e.g. 
Choud huri & Dikpatll999l) . 

Surface flux transport models (e . g. [D eVore et al] [l 984; 
Wane et al. 1989; van Ballegooiien et al. 1998; Schriiver & Title 
200 it IMackav & Lockwood l2002t iMackav et all l2004t 
Baumann et al. 2004) follow the evolution of the radial mag- 
netic field component on the stellar surface under the combined 
effects of magnetic flux emergence, differential rotation, merid- 
ional flow, and supergranular diffusion. Using empirical properties 
of bi-polar regions (i.e. latitudinal migration of emergence rates 
and tilt angles) they successfully reproduce major features of 
the solar cycle like the reversal of the polar field. Considering 
moderately rotating stars, Schriiver & Title 1 2001) showed that an 
enhancement of the rate of flux emergence produces dark polar 
caps of a single magnetic polarity , surrounded by a flu x ring of 
opposite polarity. More recently, IMackav et alJ 120041) showed 
that additional enhancements of both the latitudinal range of 
flux emergence and the meridional flow velocity are required to 
successfully reproduce a significant mixture of magnetic polarities 
at high latitudes on rapid rotators. 

The surface flux transport models may accurately describe the 
surface evolution of the magnetic field, but they do not address the 
question of whether the assumed properties of erupting magnetic 
flux are consistent with the requirement of high meridional flow 
velocities. In particular, how the latitudinal distributions of flux 
eruption are affected by enhanced rates of meridional flows. We 
therefore link our studies on the pre-eruptive and the post-eruptive 



evolution of magnetic flux to consistently quantify the impact of 
meridional circulations on the butterfly diagram and the polar mag- 
netic field properties of rapidly rotating stars. The main aim of this 
investigation is to validate the assumptions of Mackav et al. 1 2004) 
about extended latitudinal ranges of magnetic flux emergence. Yet 
the relation between observable activity properties and sub-surface 
transport mechanisms also makes it possible to infer empirical con- 
straints for specific dynamo properties. In Sect. [2] we investigate 
the rise of magnetic flux tubes through the convection zone prior 
to their emergence on the stellar surface to determine how the lati- 
tudinal eruption pattern depends on the strength of the circulation. 
The results are used in Sect. [3] to determine the impact of merid- 
ional circulations on stellar butterfly diagrams. In Sect.|4| we use 
the specified flux emergence latitudes within surface flux transport 
simulations to determine the requirements for the reproduction of 
observed stellar magnetic field properties. In Sect. [5] we discuss 
our results and their implications for possible dynamo scenarios. 
Our conclusions are summarised in Sect. [6] 



2 SUB-SURFACE EVOLUTION AND ERUPTION OF 
MAGNETIC FLUX TUBES 

2.1 Basic scenario 

Magnetic activity signatures in the atmosphere of cool stars are as- 
cribed to the eruption of magnetic flux tubes, which are generated 
by s ub-surface dynamo processes inside the convective envelope 
(e.g.[Ossendriivei 2003; Schiissler 2005, and references therein). 
The amplification of the magnetic field is expected to take place at 
the tachocline, a region of strong shear flows at the interface to the 
radiative core of the star. Helioseismological observations indicate 
that the tachocline is slightly prolate (e.g. Charbonneau et al. 1999; 
iBasu & Antiall200ll) . At the equator, the bulk of the tachocline is 
right beneath the convection zone, whereas at higher latitudes a 
substantial part of it is located inside the convection zone. In the 
dynamically unstable stratification of the convection zone, mag- 
netic flux is subject to magnetic buo yancy, w hich leads to its rapid 
loss through eruption ( Parken ll975l ; iMoreno-InsertisU l 983 ) . In the 
case of the Sun, for example, a magnetic flux tube can traverse 
the convection zone within several weeks. The generation of high 
field strengths requires magnetic fields to persist in the amplify- 
ing region over time scales comparable with those of supposed 
dynamo processes. This requirement leads to the conjecture of 
magnetic flux being stored insid e the subadiabatic ov ershoot re- 
gion beneath the convection zone l Spruit & van Ballegooiien 1982: 
van Ballegooiien 1982). Inside this stably stratified region, mag- 
netic flux tubes are perturbed through overshooting gas plumes 
penetrating from the convective envelope above. If its magnetic 
field strength is sufficiently large, a di splaced flux tube i s liable to a 
buoyancy-driven instability (cf.|Parker 1966; lD'Silva & Choudhuri 
1 1993t I Schiissler et al. 1994). Once an unstable, growing flux loop 
is properly located inside the superadiabatic convection zone, it 
rapidly rises to the stellar surface. Coriolis forces, induced by the 
internal plasma flow along the flux tube, cause an asymmetric 
evolution of the proceeding- and following legs of the tube (rel- 
ative to the direction of stellar rotation), which twists the upper 
part of the rising loop (e.g. Caligari et al. 1995). Upon eruption on 
the stellar photosphere, this twists entails a tilt between the polar- 
ity centres of the emerging bipole and the parallels on the stel- 
lar surface. Depending on their size and average field strength, 
the erupted flux tubes produce a spectrum of magnetic activity 
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signatures, from small magnetic knot s and pores to bi-po lar spot 
groups and active regions (e.g. Schri iver & ZwaarJl2000T) . Newly 
emerged surface magnetic flux dynamically disconnects from the 
sub-surface parts of the parent flux tube l Schriiver & Title 1999; 
Schiissler & Rempel 2005); in the case of the Sun this process is 
expected to take place after a few days in a depth less than 10 Mm 
below the surface. The associated upflow of entropy-rich plasma 
inside the decapitated flux tube, from its anchor rooted in the over- 
shoot region into the upper convection zone, weakens the magnetic 
field of the tube stumps and promotes their disintegration through 
turbulent diffusion driven by magneto-convective motions. The dis- 
solving magnetic field may be transported downwards to the sta- 
ble overshoot region through meridional circulation and convective 
pumping (e.g. Tobias et al. 2001), providing the seed field for fur- 
ther amplification or a new activity cycle. 

There is currently no detailed model for the removal of 'old' 
flux from the lower convection zone prior to the production of mag- 
netic flux of the opposite polarity in the successive cycle. Neither 
it is clear how, if at all, the transformation of toroidal to poloidal 
magnetic fields through the twisting of erupting flux loops is related 
to the re-generation of magnetic flux at the bottom of the convec- 
tion zone. The lack of models prevents a prediction of the amount 
and location of magnetic flux as a function of the meridional flow 
pattern and -velocities. Focusing on the transport of magnetic flux 
only, we hence presume the existence of appropriate magnetic flux 
tubes at the bottom of the convection zone of rapidly rotating stars, 
without further assumptions about their amplification and cyclic re- 
generation. 



2.2 Model setup 

2.2.1 Thin magnetic flux tubes 

The first part of the investigation is carried out in the framework of 
the thin flux tube approximation l Spruit 1981). For this approxima- 
tion to be applicable, the radius of the flux tube must be smaller 
than all other relevant length scales like, for example, the scale 
height of the pressure or the superadiabaticity, the local radius of 
curvature of the magnetic field, or the wavelength of perturbations 
propagating along the tube. Assuming an ideal plasma with infinite 
conductivity, the magnetic flux over a tube's circular cross section 
is conserved. 

Due to the relatively short signal timescale across the tube's 
diameter, the magnetic flux tube is in instantaneous pressure equi- 
librium with its environment, with the sum of the gas and magnetic 
pressure inside the tube balancing the gas pressure of the field-free 
external plasma, 
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8tt 



(1) 



The total pressure changes smoothly across the tube's surface (e.g. 
Stix 1989). For magnetic fields inside the convection zone the 
magnetic pressure is typically smaller than the gas pressure (i.e. 
B 2 /(8tt) <C p). Since the flux tube is impenetrable for the exter- 
nal plasma, perpendicular motions relative to the tube's axis cause 
the external plasma to flow around the tube. The reaction of this 
distortion on the flux tube dynamics is quantified in terms of a hy- 
drodynamic drag, 



Id 



Cd i , -, 

: Pe \Vj_\ V± 



(2) 



v±e± — (v e —v)±, between the external and internal motions, v e 
and v, respectively. The value of the dimensionless drag coefficient, 
Cd, depends on the radius, a, of the tube considered in relation 
to the spectrum of length scales of the convective motions in its 
vicinity. The value of the effective turbulent viscosity depends on 
the specific model description of the convective energy transport. 
Since it is typically much larger than the molecular viscosity, the 
drag coefficient is anticipated to be of order one. 

In addition to the interaction with its environment, realised 
through the lateral pressure equilibrium (i.e. magnetic buoyancy) 
and the hydrodynamic drag, the evolution of a magnetic flux tube 
also depends on the magnetic curvature force and the Coriolis 
force. In rapidly rotating stars, the interplay between the last two 
forces typically causes a deflection of rising flux loops to higher lat- 
itudes (e.gJChoudhuri & Gilmanll98llSchussler & Solankll992t 
ICaligari et alll995h . 



2.2.2 Stellar stratification and meridional circulation 

We consider the evolution of magnetic flux tubes in a 1 Mq 
star, described by a spherically symmetric model of the current 
Sun (i.e. R* — Rq — 6.96 ■ 10 10 cm). The outer convec- 
tion zone extends down to about 0.72 Rq and comprises at its 
lower boundary a superadiabatically stratified overshoot region of 
about 10 4 km depth. The stellar rotation period is taken to be 6 d 
(Q. = 1.212 • 10 s _1 ); for the sake of simplicity differential ro- 
tation is neglected. 

The meridional flow pattern, v p , is superposed on the stellar 
structure, anticipating that its influence on the underlying hydro- 
static stratification is negligibly small. We follow the approach of 
van Balleaooijen & Choudhuri ( 1988) and describe the meridional 
circulation analytically through the poloidal velocity components 

Vp, r = u {) (1 + Z) 2 F sin m 6 [(m + 2) cos 2 - sin 2 0] (3) 
G 
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with p e being the density of the external medium. The drag 
force tries to reduce the perpendicular velocity difference, v± = 



where F(£) and G{9) are univariate functions of the radius, r, and 
the co-latitude, 9, respectively, with £ = R*/r — 1; a more de- 
tailed description of this flow model is given in Appdx. [3] The 
circulation is parametrised through the location of its lower bound- 
ary, rt, and the dimensionless quantities n, m, and k, which also 
define the coefficients c\ and C2 in Eqs. <A5> and <A6t , respec- 
tively. We locate the maximum of the latitudinal flow velocity at 
co-latitude 6m = 53° to obtain a solar-like surface flow. Since the 
radial profiles of stellar meridional circulations are as yet virtually 
unknown, we adopt for the radial para meters rt and k the values 
used by van Balleeooiien & Choudhuri 1 1988). Thus, the assumed 
meridional flow pattern is quantified through the set of parame- 
ters ri, = 0.7 R*, n = 1.5, k — 0.5, and m = 0.76. It is com- 
posed of a single-cell circulation (per hemisphere) with a poleward 
flow at the stellar surface and an equatorward flow at the bottom 
of the convection zone (Fig. |3J. The amplitude of the circulation, 
Uo ~ —2.47 vm [cf. Eq. IA11H . depends on the choice of the peak 
flow velocity, vm, on the stellar surface. 



2.2.3 Flux tube equilibria 

The magnetic flux tubes start their evolution from a consistently de- 
fined initial configuration. For this we assume a mechanical equi- 
librium, which implies that their orientation is parallel to the equa- 
torial plane and that their radius of curvature, _Ro = ?*o cos Aq, is 
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Figure 1. Density contrast of a magnetic flux ring in mechanical equilib- 
rium at latitude Ao and depth rrj = 5.07 ■ 10 10 cm. The peak velocity of 
the meridional flow is vm = 0, 25, 50, 75, 100, 125 m ■ s — 1 (top to bot- 
tom) and the radius of the flux tube ao = 10 7 cm. 



Figure 2. Flow velocity inside a magnetic flux ring in mechanical equilib- 
rium at latitude Ao and depth ro = 5.07 • 10 10 cm. The meridional flow 
velocity is vm — 0, 25, 50, 75, 100, 125 m ■ s — 1 (top to bottom), the ra- 
dius of the flux tube ao = 10 7 cm, and the field strength Bo = 15 • 10 4 G. 
Positive (negative) values correspond to flow velocities faster (slower) than 
the stellar rotation. 



constant, where ro and Ao are the equilibrium radius and latitude, 
respectively. In the absence of meridional flows, flux rings in me- 
chanical equilibrium are non-buoyant and the magnetic curvature 
force (pointing toward the axis of stellar rotation) is balanced by 
inertia and Coriolis forces (pointing away from the rotation axis). 
The Coriolis force is caused by a plasma flow inside the flux tube 
prograde to the stellar rotation. The properties of the m echani- 
cal equilibrium are descr ibed by Spruit & van Bal legooiien ( 1982); 



Ferriz-Mas & Schiissler 1 1993, 1995); for different approache 
Choudhuri 1 1989); Moreno-Insertis et al. 1 1992l) : lFan et alj ji994^: 
Calisari et alJ<1998t) : lRempelH2003l) . 

In the presence of an equatorward meridional flow (with the 
relative direction e± perpendicular to the tube axis) the drag force, 
Eq. J2j, tries to push the flux ring to lower latitudes. For a stationary 
mechanical equilibrium to exist (Appdx.lBl. the component of the 
drag parallel to the rotation axis, e z , must be balanced by buoyancy, 
which requires the density contrast 
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(5) 



between the external and internal density, where (j e ff is the ef- 
fective gravitational acceleration defined in Eq. <B2> . The compo- 
nent of the drag force perpendicular to the rotation axis opposes 
the magnetic curvature force so that the internal flow velocity, Vo, 
required to balance the tension force, is 
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where ca is the Alfven velocity and cu the modification caused by 
the drag force in the presence of meridional circulation, given in 
Eq. <B8> . The Alfven velocity introduces the magnetic field depen- 
dence in the equilibrium condition. 

At the assumed equilibrium radius, ro = 5.07 • 10 cm, in 
the middle of the overshoot region, the stellar stratification is char- 
acterised through the pressure scale height H p = 5.52 • 10 9 cm, 
the gravitational acceleration g = 5.06 ■ 10 4 cm ■ s -2 , the density 
p e = 0.154 g ■ cm~ 3 , and the superadiabaticity 5 — -9.77- 10 -7 . 
For the meridional circulation defined above, the density contrast 
and internal flow velocity required for a stationary mechanical equi- 
librium are shown in Figs.Qand|2| respectively. 



2.3 Flux tube simulations 

The stability properties of magnetic flux tubes depend on the 
stellar stratification in its environment (including the merid- 
ional flow topology) and the magnetic field strength. The 
determination of instability criteria and characteristic growth 
times requires a specific linea r stability analysis similar to 
Spruit & van Ballegooiien ( 1982); van Ballegooijen & Choudhuri 
( 1988); Ferriz-Mas & Schiissler 1 1995), which is however beyond 
the scope of the present paper. Investigations in the case of 
the Sun and other cool stars consistently show that the evolu- 
tion of flux tubes with initial radii larger than ~ 1000 km are 
only marginally affected by the drag force, and that magnetic 
field strengths > 10 G are required to initiate rising flux loops 
whose properties upo n eruption at the surface are in agree ment 
with ob servation s (e.g . lChoudhuri & Gilmanll987tlFan et alJl994t 
Schiissler et al. 1994). Here, we consider tubes with initial radii 
ao = 70 — 1000 km and initial magnetic field strengths Bo = 
10 — 20 ■ 10 4 G, to determine the basic influence of these parame- 
ters on the latitudinal eruption pattern in the presence of meridional 
flows. The amount of transported magnetic flux, $ = Bo7rao ~ 
10 19 — 6 ■ 10 21 Mx, is compara ble to that of sm all pores and bi- 
polar spot groups, respectively iSchriiver & Zwaan 2000). We ac- 
complish simulations for the meridional flow pattern defined above 
with peak flow velocities vm < 125 m ■ s _1 . 

The simulations start with the perturbation of an equilibrium 
flux ring located in the middle of the overshoot region within the 
range of latitudes Ao = 5 — 75°. The initial perturbation consists 
of a superposition of harmonic displacements with low wave num- 
bers and amplitudes of a few percent of the local pressure scale 
height. The adiabatic evolution of individual flux tubes is followed 
using a non-linear Lagrangian scheme described in more detail in 
More no-Insertisl I 1986) and Cali gari et alHl995h . Owing to the lat- 
eral pressure balance, Eq. Q, the summit of a rising flux loop ex- 
pands as the external pressure decreases. The simulations stop just 
beneath the stellar surface (r > 0.98 R*), where the underlying thin 
flux tube approximation becomes inapplicable. In the upper-most 
part of the convection zone the rise of the tube summit is domi- 
nated by strong magnetic buoyancy so that the tube trajectory is 
almost radial. The final point of the simulations thus determines 
the location of the emergence of a bi-polar spot group. 
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Figure 3. Latitudinal trajectories of the summit of rising flux loops with 
initial field strength Bo = 15 • 10 4 G and tube radius ao = 10 7 cm 
for maximal flow velocities vm = 10, 75, and 100 m ■ s — 1 (grey shaded 
lines). The flow of the meridional circulation (solid lines, n = 1.5, m = 
0.76, k = 0.5, Tf, = 0.7 R*) is poleward at the surface and equatorward at 
the bottom of the convection zone. The shaded region at the bottom of the 
convection zone marks the overshoot region. 



2.4 Latitudinal distributions of magnetic flux eruption 

The interplay between Coriolis and magnetic tension force results 
in a deflection of rising flux loops to higher latitudes. Without 
a meridional circulation the deflection is less than AA = A e — 
Ao < 10°, with A e being the eruption latitude on the stellar sur- 
face. The deflection is greatest for flux tubes which originate from 
low initial latitudes and decreases for higher starting latitudes. In 
the presence of meridional flows the poleward deflection is consid- 
erably larger and shows a characteristic dependence on the initial 
latitude of the flux ring (Fig.|5J. Flux tubes originating from low to 
intermediate latitudes (here, around Ao ~ 25°) are most affected 
by poloidal flows, whereas those starting from equatorial or high 
latitudes are in contrast only slightly more deflected than in the ab- 
sence of an external circulation (Fig. |4}- This latitude-dependent 
deflection resembles closely the initial density contrast and internal 
velocity variation (Figs.Qand|2| respectively), suggesting that the 
strong deflection is partly caused by the equilibrium conditions. 

The latitudinal distribution patterns obtained for different 
meridional flow velocities are qualitatively similar, with the ampli- 
tude of the pattern increasing with the peak flow velocity (Figs.|4^ 
and[5). Whereas flux tubes with smaller radii show a characteristic 
deflection pattern with a maximum at intermediate latitudes, flux 
tubes with radii larger than 200 — 300 km yield an eruption pat- 
tern in which the poleward deflection continuously decreases from 
low to high initial latitudes (Fig.|4j)). Flux tubes with larger initial 
radii show deflections only slightly larger than in the case without 
a meridional circulation. For magnetic flux tubes with high field 
strengths the latitudinal deflection is weaker (Fig. |4j;), since the 
relative impact of the magnetic buoyancy on the dynamical evolu- 
tion is larger. Their trajectory through the convection zone is more 
radial, with minor latitudinal deflections only. 




Figure 4. Latitudinal deflection, A A = A e — Ao, of erupting flux tubes in 
the presence of meridional flows with different peak flow velocities (Panel 
a, for Bo = 15 ■ 10 4 G and ao = 10 7 cm); initial tube radii (Panel b, for 
vm = 100 m ■ s" 1 and Bo = 15 ■ 10 4 G); initial magnetic field strengths 
(Panel b, for vm = 100 m ■ s — 1 and ao = 10 7 cm). Gaps in the curves 
for high flow velocities and low magnetic field strengths are due to stable 
flux ring equilibria for the respective set of parameters. 
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Figure 5. Dependence of the poleward deflection, AA = A e — Ao, on the 
meridional flow velocity, Vm ■ 



2.5 Comments 

Thin magnetic flux tubes with weak average field strengths are 
found to be susceptible to meridional circulations. In particular 
flux tubes originating from low to intermediate latitudes are sub- 
ject to significant poleward deflections during their rise to the stellar 
surface. We conjecture that this latitude-, size-, and field strength- 
depending deflection causes the bulk of low-flux elements to erupt 
at higher latitudes than more substantial flux tubes. Compared to 
stars with weak or no meridional flows, this separation effect causes 
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Table 1. Correspondence between the latitudes of origin, Ao, and latitudes 
of emergence, A e , of flux tubes with an initial field strength of Bq = 15 ■ 
10 4 G and radius ao = 100 km (see Fig. 161. 
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-10 
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-9 


60 


-9 


L 


58- 
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- 10 
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-9 


63 


-9 


65 


-9 



a bias in the overall surface distribution of erupting magnetic flux 
toward higher latitudes. 

The enhanced susceptibility of smaller magnetic flux elements 
inside the convection zone to meridional flows is in basic agree- 
ment with the properties of magnetic flux features observed on 
the solar surface. There, smaller flux elements are readily trans- 
ported to higher latitudes, whereas new active regions and sunspots 
(in their entity) hardly participate in the poleward motion (vm ~ 
10 m ■ s^ 1 ). However, after their fragmentation through supergran- 
ular motions, the dissolving remnants of initially large flux concen- 
trations are swept polewards as well. 



3 STELLAR BUTTERFLY DIAGRAMS 

The latitudinal distribution of erupting flux tubes depends in a 
non-linear way on the magnitude of the meridional circulation as 
well as on the magnetic flux and original latitude of the tube in- 
side the overshoot region. In the following, we consider the emer- 
gence pattern of flux tubes with an initial magnetic field strength 
Bo = 15 ■ 10 4 G and an initial radius ao = 100 km as representa- 
tive for the full range of parameters considered in Sect.|2| 

Surface flux transport models require a description of the sta- 
tistical properties of newly emerging bipoles during an activity cy- 
cle, including their emergence rates, latitudes, sizes, fluxes, and tilt 
angles. Since the observational data base on rapidly rotating solar- 
like stars is insufficient to create empirical butterfly diagrams, we 
adopt the solar butterfly diagram as a qualitative template, which 
we extrapolate in latitude and scale according to the results in Sect. 
|2| The solar butterfly diagram is characterised by emergence lati- 
tudes between about 40° latitude at the start of a cycle and 5° at 
its end. The width of the wing of each cycle is about 10° at the 
beginning of the cycle and decreases to 5° at the end. The cycle 

i ieriod is 11 yr, e xcluding one-year overlaps of successive cycles 
van Ballegooiien et all l998: Mac kay et all2004) . 



3.1 Latitudes of emergence vs. latitudes of origin 

Based on the simulation results in Sect. [2] we associate latitudi- 
nal ranges of magnetic flux emergence, A e , with latitudinal ranges 
of origin, Ao, of the parent flux tubes in the overshoot region. For 
the reference (i.e. solar) butterfly diagram, we refer to the case 
Vm = 0m • s _1 , since the weak solar meridional circulation with 
a peak velocity of ~ lOm-s -1 has little effect on the sub-surface 
evolution of rising flux tubes. In this case, the relationship between 
A e and Ao is virtually linear (Fig. [6^, solid line/asterisks), and 
we ascribe the emergence of bipolar regions on the Sun to mag- 
netic flux tubes which originate from within the range of latitudes 
Ao — 35 — 5° . Further assignments between latitudes of emergence 
and latitudes of origin are summarised in TableQ We shall refer to 
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Figure 6. Panel a: Latitudes of emergence, A e , of magnetic flux tubes 
as a function of their latitude of origin, Ao, inside the convection zone. 
The initial magnetic field strength of the flux tubes is Bq = 15 • 10 4 G 
and the initial radius ao = 100 km. The flow velocity of the merid- 
ional circulation is Vm = Oms^ 1 (solid), 50m ■ s _1 (long-dashed), 
75 m ■ s _1 (short-dashed), 100 m • s -1 (dashed-dotted), and 110 m • s _1 
(dashed-tripple dotted); lines represent polynomial fits to the results (sym- 
bols) of the flux eruption simulations in Sect. |2| Panel b: Deflection ra- 
tio, A e (^M )/A e (0), relative to the case without any meridional circulation. 
Values larger (smaller) than one indicate relative deflections toward the pole 
(equator). Note that the abscissa in panel b is equal to the ordinate in panel 
a, since the deflection is expressed in terms of latitudes of emergence. 



the small (Ao = 35-5°), intermediate (47-3°), and large (58-3° 
ranges of originating latitudes as case S, I, and L, respectively. 



3.2 The wings of the butterfly 

In contrast to the solar reference case, strong meridional circula- 
tions imply non-linear relationships between the eruption latitudes 
of magnetic flux tubes and the latitudes of their origin (Fig.[3Ji) and, 
consequently, distortions of the wing shapes of stellar butterfly di- 
agrams. To transform the solar template into a stellar butterfly di- 
agram, we determine the ratio between the latitudes of emergence, 
Xb(vm), subject to meridional flows and the latitudes of emer- 
gence, A e (0), of the reference case with vanishing meridional flow. 
The deflection ratios A e (vM)/Ae(0) in Fig.[6j) represent mappings, 
which describe the non-linear stretching of the solar template. In 
addition to the poleward displacement of flux emergence at inter- 
mediate latitudes also a slight equatorward displacement occurs at 
low latitudes. 

Examples of stellar butterfly diagrams subject to meridional 
circulations with different peak velocities are shown in Fig. Q 
(for case I). Figure^ depicts a stellar butterfly diagram subject 
to a weak meridional circulation (i.e. the solar template linearly 
stretched to the latitudinal range 50 — 10°), whereas the Figs.0) & 
c show the result of the non-linear transformation of the latitudinal 
pattern of emerging bipoles caused by higher flow velocities. The 
diagrams show that for strong meridional circulations the number 
of bipoles emerging at higher latitudes is increased: the higher the 
flow velocity, the more convex is the temporal evolution of the stel- 
lar butterfly diagram. Qualitatively similar results are obtained for 
the cases S and L. 
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Figure 7. Butterfly diagrams over four 1 1 -year activity cycles in the pres- 
ence of meridional circulations with flow velocities vm = 10 m ■ s _1 
(Panel a); vm = 75 m ■ s _1 (Panel b); Vm = 100 m ■ s _1 (Panel c). 
In all cases the generating magnetic flux tubes originate from within the 
range of latitudes Ao = 47 — 3°. The pre-emptive poleward deflection 
caused by the meridional circulation pushes the upper boundary of the but- 
terfly diagram to higher latitudes, from A e » 50° for the weak flow (Panel 
a) to A e « 60° for the strong flow (Panel c). 



of flux with the same (opposite) polarity through the convective 
turnover of supergranular cells leads to an enhancement (annihi- 
lation) of magnetic flux. Averaged over length scales larger than 
supergranular cells (in the case of the Sun about 30 Mm) the evolu- 
tionof the magnetic field resembles a dispersion process I Leiehton 
1964). The lifetime of individual magnetic surface features depends 
on the diffusion time scale. The evolution of the global magnetic 
field topology, in turn, depends on the relation between the lifetimes 
of merging flux features and the characteristic transport timescales 
of the differential rotation and the meridional circulation. Differen- 
tial rotation shears and diverges bipolar spot groups spread out over 
a range of latitudes. The increasing distance between the centres of 
opposite flux polarities reduces the rate of mutual flux annihilation, 
whereas the associated decrease of the tilt angle reduces the bias in 
the latitudinal distribution of magnetic polarities (see also Appdx. 
O . Merging magnetic surface features are transported to higher lat- 
itudes by the poleward directed meridional circulation and pile-up 
in polar regions. Owing to the biased latitudinal polarity distribu- 
tion, the polar flux is dominated by one polarity during half of the 
magnetic cycle. With the reversal of the polarity of newly emerging 
bipolar groups after half a magnetic cycle (Hale's law), magnetic 
flux of the opposite polarity is predominantly transported toward 
the pole, which annihilates and substitutes the previous flux con- 
glomeration. This alternating process implies a stable oscillation of 
the net surface flux over successive activity cycles, without a per- 
sistent accumulation of magnetic flux. For a review about the devel- 
opment and limitations of surface flux transport model see Sheelev 
(e.g.Hxi). 

In contrast to the case of the Sun, where the polar mag- 
netic field is virtually unipolar, with magnetic flux densities of 
about 10Mx-cm -2 , the surface magnetic fields observed on 
young active stars exhibit a high degree of intermingling of mag- 
netic fluxpolarities toonati & Co llier Cameron 1997; Dona ti et alj 
1999, 2003) and dark polar caps (Strassmeier 2002), indicative for 
strong magnetic field strengths. Schriiver & Title 1 2001 ) found that 
flux emergence rates thirty times larger than on the Sun entail polar 
flux densities of 300—500 Mx ■ cm" 2 , which are deemed to be suf- 
ficiently large to suppress convective upflows of energy, entailing 
the formation of dark spots (e.g. Priest 1982). In their simulations, 
the polar magnetic flux is unipolar and surrounded at lower lati- 
tudes by a ring of magnetic flux of the opposite polar ity. Focusin 
on the intermingling of polar magnetic flux, JMackav et all 1 200' 
identified two key parameters for the generation of high degrees of 
polarity mixture: magnetic bipoles have to emerge at latitudes up to 
50 — 70°, and the peak value of the meridional flow has to be over 
100 m-s"" 



The simulations of Mackav et alj l2004h disregard the feed- 
back of strong meridional circulations on the latitudinal pattern of 
flux eruption. This feedback is now provided through the consis- 
tently determined s tellar but terfly diagrams in Sect. [3] In the fol- 
lowing, we use the Macka v et alJ -model to investigate the impact 
of the convex wing structures on the polar magnetic field proper- 
ties. 



4 SURFACE EVOLUTION AND INTERMINGLING OF 
MAGNETIC BIPOLES 

4.1 Basic scenario and empirical constraints 

After the dynamical disconnection of newly emerged bipoles from 
their parent flux tubes, the small-scale evolution of magnetic sur- 
face features is governed by local magneto-convective motions and 
interactions with ambient magnetic flux concentrations. Advection 



4.2 Surface flux transport model 

The surface flux transport simulations of Ma ckav et alj 120041) 
follow the radial magnetic field component on the stellar sur- 
face regarding the combined influence of flux emergence, dif- 
ferential rotation, meridional flow, and supergranular diffusion. 
Using spherical harmonics up to degree 63, the spatial res- 
olution is about 30 Mm. The simulations are based on the 
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extrapolatio n of solar transport and flux emerging prop erties 
jGaizauskas et alJll983tlWang & Sheelevlll989t [Harvey & Zwaanl 
ll993USchriiver & Harvevlll994lTian et allll999f> to more rapidly 
rotating stars (see also Schriiver & Title 2001). 

The strength and profile of the differential rotation of rapidly 
rotati ng stars is similar to the solar case (e.g. Coll ier Cameron et alj 
2002). We therefore adop t the solar differential rotation profile 
given by Snodgrass ( 1983), which implies the characteristic shear 
timescale tq — 0.25 yr. The meridional flow profile is given by 
fp(A) = —VMsm(n\/X v ), whereb y the flow v elocity vanishes 
near the pole above X p = 75° lHathawav 1996). The character- 
istic timescale for the latitudinal flux transport is r m f = Rq/vm- 
Lacking appropriate empirical constraints, it is assumed that the su- 
pergranular convective profile and turn-over timescales of rapidly 
rotating stars are so lar-like, implying a di ffusion coefficient of 
D — 450 km 2 ■ s _1 J van Ba llegooiien et al. 1998) and a diffusion 
time scale for magnetic features with length scale I of tb = I 2 /D 
(= 34 yr for I = Jfo). 

Following Schriiver & Title (2001), we take the net amount of 
magnetic flux emerging during each cycle to be thirty times the so- 
lar value, that is 3.15 • 10 26 Mx. New flux is injected into the evolv- 
ing surface magnetic field distribution at random longitudes in the 
form of 13 200 bipolar mag netic regi ons, whose average tilt angles 
vary as A e /2 l Wang & Sheelev 1989). In each successive cycle the 
polarity of the preceding and following flux of the bipoles alter- 
nates according to Hale's law. The simulations cover a time span 
of four activity cycles to verify that stable oscillations of the po- 
lar field from one cycle to the next are obtained. For further details 
about model assumptions, input parameters, and initial conditions 
see lMackav etail<2004 Sect. 2-4). 

In the following, we focus on the total (unsigned) magnetic 
flux, 



$tot(t)= / \B r (R,0, 4>,t)\dS, (7) 

the total (unsigned) polar magnetic flux in each hemisphere, 

*[np,s P ](t) =J \B r (R,e [g<20 o e>160 o y <p,t)\dS, (8) 

the imbalance of positive and negative magnetic flux in the polar 
regions, 



w*)=y i ww,*> 1 6o°]^.') d5 ' ( 9 ) 

and the polar magnetic flux in the northern hemisphere, 
* [+H (t)= / ' B T (R,e <20 o ,<l>,t) [Br>0 , Br<0] dS • (1Q) 



The principal quantities of the analysis will be mean magnetic flux 
densities, that is the respective magnetic flux divided by the sur- 
face over which it is calculated (see Mackav et al. 2004). The last 
quantity, Eq. ilOt is a measure for the degree of intermingling of 
opposite polarity elements at high latitudes. 



4.3 Surface magnetic field properties of rapid rotators 

We consider weak and strong meridional flows of the case I to de- 
tail our results (Fig. [8). For both flow velocities the mean surface 
flux density, <I>tot/(47r7? 2 ), is in phase with the activity cycle emer- 
gence rate (Fig. [8^ & f). If the meridional flow velocity is high, 
then more magnetic flux is pushed to higher latitudes, where the 
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Figure 10. Panel a: Cycle-averaged ratio between the positive and negative 
magnetic flux within the northern polar cap as a function of the meridional 
flow velocity, vm, for originating latitudes Ao = 35 — 5° (triangles); 47 — 
3° (diamonds); and 57 — 3° (asterisks). Panel b: Cycle-averaged polar flux 
densities. 



proximity of opposite bipole polarities increases the flux annihi- 
lation rate. In this case, the mean surface flux densities are conse- 
quently smaller than for slow meridional circulations. The resulting 
variation of the mean polar flux density (i.e. beyond 70° latitude) 
is shown in Figs. [8p & g. In the case of weak meridional flows 
this quantity shows a noticeable phase shift relative to the mean 
surface flux density, with a cycle-averaged polar flux density of 
~ 192 G. In the case of fast flows the polar field rises and falls in 
phase, with a significantly higher cycle-averaged polar flux density 
of ~ 270 G. Breaking up the polar flux density imbalance (Figs. 
[SJ; & h) according to the contributions of positive and negative flux 
(Fig.|8ji & i) shows that for slow meridional flow velocities the po- 
lar field is unipolar, with a field reversal midway through the cycle. 
For fast flow velocities, in contrast, both magnetic polarities are 
present throughout the cycle and rise and fall in strength together. 
The ratio between the positive and negative flux densities in the 
polar caps is thus in both cases different (Fig. [8^ & j). For weak 
flows there is hardly any intermingling of flux of opposite polar- 
ities at high latitudes, except for a limited time span during field 
reversal, when little magnetic flux is located within the polar cap. 
For strong meridional flows, in contrast, the intermingling of mag- 
netic flux with opposite polarities is about 65%, averaged over an 
activity cycle. 

Snapshots of maps showing the radial magnetic field compo- 
nent (Fig. [5) illustrate that in the case of weak flows the magnetic 
flux in polar regions is dominated by one polarity only, whereas in 
the case of fast flows the field is intermingled. 



4. 3. 1 Parameter study 

For the latitudinal ranges I and L, the degree of intermingling of 
opposite polarities within the polar regions exceeds 20% for merid- 
ional flow velocities beyond about 40 m ■ s _1 (Fig. llOfc ). The fact 
that in case I the intermingling is higher than in case L indicates 
that merely increasing both the meridional flow velocity and the 
latitudinal range of flux origin does not produce per se higher de- 
grees of intermingling. Instead, there exists an optimal combination 
of latitudes of origin and meridional flow velocities for which the 
intermingling within polar caps peaks. 
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Figure 8. Results of the surface flux transport simulations for meridional circulations with a weak flow velocity of 11 m • s — 1 (Panels a-e), and a strong flow 
velocity of 100 m ■ s _1 (Panels f-j); the associated butterfly diagrams are shown in Figs.0 and0;, respectively, with generating flux tubes originating from 
within the range of latitudes Ao =47 — 3° (Case I). Panels a & e: mean surface flux density, <J>t ot /(47ri? 2 ); panels b & g: mean polar flux density of the 
northern (solid) and southern (dotted) hemispheres; panels c & h: signed polar flux densities in the northern (solid) and southern (dotted) hemisphere, with 
dashed lines representing the initial values; panels d & i: positive (solid) and negative (doited) polar flux density of the northern hemisphere; panels e & j: 
ratio between the positive and negative polar flux density of the northern hemisphere. Each quantity is determined once every 27 days. Dashed-dotted lines 
mark the start of a new cycle. All flux densities are in Mx ■ cm -2 . 



The associated polar flux densities are on the order of 200 — 
300 Mx • cm~ 2 (Fig. 110b ). with values increasing (decreasing) 
with the meridional flow velocity above (below) vm ~ 50 m • s _1 . 
This behaviour is due to the different transport- and annihilation- 
timescales (see Appdx.lcl. In the case of strong meridional flows, 
the advection of magnetic flux to higher latitudes is more efficient 
than the decrease caused by flux annihilation as the magnetic sur- 
face features are transported poleward. The same argument holds 
for the increase of the polar flux density with larger latitudinal 
ranges of flux origin, since for bipoles emerging at higher latitudes 
the shorter migration time towards the pole implies less flux annihi- 
lation. For meridional flow velocities around 100 m ■ s _1 the polar 
flux densities become sufficiently large (~ 300 G) to suppress con- 
vection. 

In case S, the polar magnetic flux is unipolar throughout the 
activity cycle, without a significant intermingling of opposite po- 
larities. The average polar flux densities of fa 300 G are almost 
independent of the meridional flow velocity. 

In summary, we find that a combination of intermediate ranges 
of originating latitudes (here, Ao =47 — 3°) with fast meridional 



flow velocities (vm ~ 100 m • s _1 ) is optimal for the production 
of both strong intermingling of polarities and high field strengths in 
polar regions. This suggests that in rapidly rotating solar-like stars 
magnetic flux originates from latitudinal ranges somewhat larger 
than in the case of the Sun. 



4.3.2 Effect of the pre-emptive poleward deflection 

Mack av et alJ J2004I) disregard the influence of meridional circu- 
lations on the sub-surface evolution of magnetic flux. In their in- 
vestigation the reference (i.e. solar) butterfly diagram is linearly 
stretched to different latitudinal ranges, leaving the structure of the 
wings unaltered. According to our results in Sect. 13.21 the eruption 
pattern within a given range of emergence latitudes depends on the 
meridional flow velocity, with more bipoles emerging at higher lat- 
itudes the faster the flow velocity. We compare the two approaches 
in the case of a meridional circulation with a flow velocity of 
100 m • s" 1 . Keeping other model parameters unchanged, includ- 
ing the pre-eruptive poleward deflection entails about 20% higher 
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Figure 9. Maps of the radial magnetic field component in the case of weak (fO m ■ s — , left column) and strong (fOO m • s — , right column) meridional 
flow velocities. The associated butterfly diagrams (cf. Fig.0 & c) are generated by magnetic flux tubes which originate from within the latitudinal range 
Aq = 47 — 3°. Each set shows snapshots after year 1 (cycle minimum, Panels a & e), 5 (before maximum, Panels b & f), 8 (after maximum, Panels c 
& g), and 1 1 (declining phase, Panels d & h). White and black shadings represent positive and negative magnetic flux densities, respectively, saturating at 
±300Mx-cm- 2 . 



polar flux densities than in the original approach of Macka v et alJ 
This result applies to all ranges of generating latitudes (Table |2j. 
The reason for this increase lies in the more convex structure of 
the wings of the butterfly (cf. Fig.[7}. The larger number of emerg- 
ing bipoles at higher latitudes over an activity cycle implies shorter 



transport times (on average) toward the pole during which less flux 
annihilation occurs. 
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Table 2. Comparison of the cycle-averaged polar magnetic field strength, 
B, including and excluding the effect of pre-eruptive deflection of mag- 
netic flux tubes for a meridional circulation with flow velocity vm = 
f 00 m ■ s~ 1 . The latter values are taken from Mackav et al. 1 2004). 
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5 DISCUSSION 

The present work relates observed properties of magnetic flux on 
the stellar surface to the sub-surface evolution of the originating 
magnetic fields and, therefore, to the underlying dynamo processes 
inside the convectiv e envelope. Our results extend the findings of 
Mackav et al. (2004), who identified high latitudes of magnetic flux 
emergence and high meridional flow velocities as the key ingre- 
dients for a strong mixing of opposite polarities in polar regions 
of rapidly rotating stars. We find that the enhanced pre-eruptive 
poleward deflection of rising flux tubes inside the convection zone 
provides a consistent explanation for the required larger latitudinal 
range for newly emerging bipoles. 

Within the framework of our model assumptions, a solar-like 
butterfly diagram with flux emergence up to 40° latitudes is gen- 
erated by erupting magnetic flux tubes, which originate at the bot- 
tom of the convection zone from latitudes < 35° . Owing to the 
slight prolateness of the solar tachocline (e.g. Charbonne au et alj 
1999), these low latitudes coincide in the case of the Sun well with 
the region of strong shear flows, efficient magnetic field amplifica- 
tion, and field storage. For more rapidly rotating stars, the butterfly 
structures producing the optimal levels of polar flux intermingling 
are generated by erupting flux tubes, which originate from higher 
latitudes, < 50°, as well. This suggest the existence of a larger lat- 
itudinal range of efficient dynamo operation than in the case of the 
Sun, implying a possibly less prolate tachocline. 

Based on the results of Schriiver & Title (2001), we assumed a 
magnetic flux emergence rate thirty times larger than the solar value 
to produce polar flux densities sufficiently strong to entail dark 
caps. Dynamo operation is generally expected to be more efficient 
at higher rotation rates, mainly entailing an increase of the surface 
filling factor of magnetic features (e,g. lSaaJll99ll) . The larger lat- 
itudinal range of dynamo operation suggested by our results com- 
plies with this expectation of a larger amount of magnetic flux gen- 
erated during an activity cycle. Considering a possible concomitant 
increase of the average magnetic field strength, we note that strong- 
field flux tubes are less susceptible to meridional circulations and 
therefore hamper the formation of high flux densities and high de- 
grees of flux intermingling in polar regions. We also find that flux 
tubes with initial radii ao > 1000 km, which are expected to lead 
to the formation of average sunspots, are only marginally affected 
by meridional flows. Since the magnitude of the drag force scales 
oc Vj_/a, a meridional flow velocity about three times faster would 
correspond roughly to a flux tube radius about ten times smaller, 
entailing a noticeable change of the eruption latitude. In view of 
the uncertainties of the sub-surface flow profile, this estimate indi- 
cates that also larger magnetic flux elements could be affected by 
meridional circulations. 

By using the solar butterfly diagram as an empirical template, 
we circumvent the treatment of the magnetic field amplification and 
dissolving mechanisms inside the stellar convection zone, which as 



yet lack a consistent theoretical description. This approach, how- 
ever, implies that the activity signatures of rapidly rotating stars 
are due to a solar-like dynamo. The assumption of an afi-dynamo 
mechanism has important implications, for example, on the period, 
and even the mere existence, of the activity cycle. The 11 yr activity 
cycle assumed here sets the time scale for the pile-up of magnetic 
flux in polar regions. If we, for example, increase the cycle length, 
then we would obtain qualitatively similar results by increasing the 
total amount of magnetic flux emerging per cycle. A conclusive pa- 
rameter study of different ratios between eruption, transport, and 
diffusion time scales is however beyond the scope of this work; a 
discussion of the impa ct of selected su rface transport model param- 
eters is given by Ma ckav et al J 12004). 

Our work supplements earlier investigations, which con- 
sidered the evolution of the weak solar magnetic field subject 
to a prescribed meridional flow pattern and a dynamo wave 
at the bottom of the convection zone (e.g Dikoati & Choudhuri 
1994; Choudhuri et al. 1995; Choudhuri & Dikp atil 19991) . All ap- 
proaches, however, suffer from the same lack of knowledge about 
the sub-surface profiles of meridional circulation and their depen- 
dence on the rotation, mass, structure, and evolutionary stage of the 
star. Owing to the density difference between the top and the bot- 
tom of the convection zone, the equatorward flow is slower than 
the surface flow, and the larger the radial extent of the equatorward 
flow region, the lower the associated velocities will be. The latitu- 
dinal flow velocity at the bottom of the convection zone possibly 
has an influence on the propagation velocity of dynamo waves and, 
consequently, on the period of the activity cycle. But in absence of 
an appropriate theory, this impact cannot be quantified. The circu- 
lation model applied here assumes a ratio of 2:1 between the up- 
per convection zone flowing poleward and the lower convection 
zone flowing toward the equator. A different radial profile would 
not only affect the equator- and poleward drift of rising flux tubes 
on their way to the surface, but also their equilibrium and stability 
properties inside the overshoot region. The former aspect has to be 
investigated in a parameter study considering both meridional flow 
and differential rotation profiles as well as different stellar rotation 
periods. Furthermore, specifically determined tilt angles and erup- 
tion time scales of newly emerging bipoles have to be consistently 
included in the surface flux transport model. 

Compared with a mere linear stretching of the solar butter- 
fly diagram to higher latitudes (e.g. Mackav et al. 2004), the addi- 
tional convexity of the butterfly wings promotes the intermingling 
of polar flux even more, which eases the necessity for high pole- 
ward flow velocities. Although the meridional flow velocities of ~ 
100 m ■ s" 1 are lower than the values suggested by Mack av et alJ 
1 2004), they are yet still about ten times larger than in the case of the 
Sun. Whereas in the case of the Sun, with vu ~ 11m • s, the 
bias in the latitudinal size distribution of newly emerging bipoles 
is probably neither discernible nor relevant, convex butterfly struc- 
tures and the size distributions of emerging dipoles may become 
observable on rapidly rotating stars once the temporal and spa- 
tial resolutions are sufficiently high. Strong meridional circulations 
on the order of 100 m • s _1 are also expected to cause a shift of 
the stellar eigenfrequencies, in addition to the splitting of eigen- 
frequencies caused by differential rotation (M Roth, priv. comm.) 
Depending on the sub-surface profile of the flow, a meridional cir- 
culation of 100 m • s" 1 would cause in the case of the Sun a fre- 
quency shift of a few tenth of /^Hz. Albeit small, this effect would 
be observable on the Sun. Given sufficiently long and complete 
time sequences, astroseismological observations may provide con- 
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straints on the magnitude (and possible structure) of stellar merid- 
ional flows. 



6 CONCLUSION 

Strong meridional circulations enhance the poleward deflection of 
rising magnetic flux tubes inside stellar convection zones. Flux 
tubes which comprise small amounts of magnetic flux and which 
originate between low to intermediate latitudes within the over- 
shoot region are particularly susceptible to meridional flows. The 
resulting latitudinal distribution of newly emerging bipoles ren- 
ders the wing structure of stellar butterfly diagrams distinctively 
more convex than in the solar case. The larger amount of magnetic 
flux emerging at higher latitudes supports the formation of high 
magnetic flux densities and a significant intermingling of opposite- 
polarity magnetic flux in polar regions, which is in agreement with 
recent observational results of rapidly rotating solar-like stars. 

The strong pre-eruptive deflection of magnetic flux tubes pro- 
vides a consistent explanation for the required high latit udes of 
flux e mergence identified by previous investigations l Mack ay et alj 
2004). The synergetic combination of pre-eruptive flux tube deflec- 
tion and post-eruptive bipole transport yields higher values for the 
average polar flux densities than in the previous approach, and thus 
eases the necessity for high meridional flow velocities. The func- 
tional dependence of the polar field properties on the extent and 
structure of the stellar butterfly diagram makes it possible to con- 
jecture about potential regions of efficient dynamo operation at the 
bottom of the convection zone. In rapidly rotating cool stars, we 
suppose the latitudinal range of magnetic flux production to extend 
to higher latitudes ( < 50°) than in the case of the Sun ( < 35°). 
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APPENDIX A: MERIDIONAL FLOW MODEL 

We u se the analytical model of van Ballesooijen & Choudhuri 
1 1988) for a parametrised description of the meridional circulatio n 
inside the convection zone (see also;Dikpati & Choudhuri 19941) . 
The poloidal flow pattern, v p = (V x \l/eV) /p e , is expressed in 
terms of the scalar function, ^/, to ensure that the stationary conti- 
nuity equation is consistently fulfilled. Under the assumption of a 
separation of variables, the stream function 

* (r,9)rsm9 = u Rlp e {r) ■ F{r) ■ G (6) , 



is expressed in terms of the univariate functions 
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depend on the location of the lower boundary, £(, = R*/rt — 1, of 
the flow pattern, with r% being a free model parameter. The poloidal 
components of the flow velocity are 
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Figure Al. Radial (left) and latitudinal (right) velocity components of the 
meridional flow for the parameter n = 1.5, m = 0.76, k = 0.5, and 
r\y = 0.7/?*. The maximum of the latitudinal flow velocity on the stellar 
surface is located at 9m = 53°. At the radius ~ 0.87?*, the latitudinal flow 
pattern changes its sign from poleward (above) to equatorward (below). 
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and the local density scale height H p = — (d In p e /dr ) . If the den- 
sity stratification is approximated through the expression p e (r) = 
Cf n , then it is TL = n and Eq. <A9> results in Eq. J4j. For an adi- 
abatically stratified convection zone with a ratio of specific heats 
7 = 5/3, it is n — 1.5 (van Ballegooijen & Choudhuri 1988). 
The maximum of the latitudinal flow velocity is located at co- 
latitude 9m, specified by the condition tan#M = V m + 1- We 
use the value m = 0.76 to obtain a solar-like value 9m — 53°. 
The flow pattern is composed of a downflow at high latitudes (i.e. 
v p ,r < for < 9 < 58°) and an upflow at low latitudes (i.e. 
V PlT > for 58° < 9 < 90°). The latitudinal velocity compo- 
nent changes its direction at about 0.8i?*, that is v P: e > (equa- 
torward) for 0.7 < r//?*<0.8 and v P: e < (poleward) for 
0.8<r/E* < 1 (Fig. ED. At the surface, r = R+ it is £ = 
and Vv.e — ~uo sin m+1 9 cos 9. The amplitude, 
Vm 



v p ,e 



Uo — — - 



(All) 



fl 9 m cos 9m 
of the circulation is determined through the maximal flow velocity, 
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APPENDIX B: FLUX TUBE EQUILIBRIUM 

The mechanical equilibrium in the presence of meridional flows is 
determined in the frame of reference co-rotating with the star. For 
a homogeneous flux tube (i.e. vanishing derivatives with respect to 
the arc length), the stationary equation of motion is 



-^ 



B- 



-i-irpj R 



B' 2 



V [p + —J +pg eff +2p (vxti) , (Bl) 



where v is the internal flow velocity, B the magnetic field strength, 
p the density, p the gas pressure, R the local radius of curvature, n 
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the normal vector of the tube, Q, = f2e z the stellar rotation vector, 
and 
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the effective gravitational acceleration comprising the centrifugal 
contribution, which (for the model assumptions specified in Sect. 
12.2.31 inside the overshoot region is of the order O (Q. 2 r/g) ~ 
10~ 4 . Anticipating the meridional circulation to be of minor im- 
portance for the stellar structure, the external stratification is ap- 
proximately hydrostatic, 
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Using the lateral pressure balance, Eq. 0, the difference between 
Eqs. OH) and EJ yields 
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with the hydrodynamic drag force, fo , accounting for the influence 
of the distorted external flow on the dynamics of the flux tube. For 
toroidal magnetic flux tubes with a constant radius of curvature, the 
component of Eq. <B4t tangential to the tube axis (i.e. the azimuthal 
component) vanishes per se. The bi-normal component (i.e. parallel 
to axis of stellar rotation, e z ) implies the condition 
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With the definition for fu in Eq. J2j this yields the density con- 
trast given in Eq. J5J, which is required to balance the influence of 
the drag force along the z-axis through buoyancy. The component 
perpendicular to the stellar rotation axis (e^ = — ft), 



50 100 150 200 

Time, [days] 
(B) 



250 300 



Anp R^\ p 



9eff ■ en) 



-2vCl 



-- (/d -en) 
p v i 



Figure CI. Panel a: North polar flux density for a single bipole inserted at 
a mean latitude of emergence, A = 21° (dashed), 39° (dotted), and 44° 
(solid). Panel b: Ratio between positive and negative magnetic flux in the 
polar region; the closer to one, the stronger the intermingling. 
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shows that a meridional flow toward the equator in the overshoot 
region supports the outward directed inertia and Coriolis force in 
balancing the magnetic tension force of the curved flux ring, where 
R = r cos A is the constant radius of curvature. From the alterna- 
tive form of Eq. I B 6i . 

^ + 2 ^ + t = f (B?) 

follows the internal flow velocity given in Eq. J6|, which is required 
to keep the flux ring in mechanical equilibrium. The contribution by 
the drag force is formally expressed in terms of the velocity 
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whereas the Alfven velocity <? A = B 2 /(4-Kp) contains the depen- 
dence on the magnetic field strength. 

The particular susceptibility of magnetic flux rings located at 
low to intermediate latitudes (cf. Figs.Q&|2j is likely a character- 
istic feature of the equilibrium properties. Neglecting centrifugal 
forces, the latitudinal variation of the density contrast, Eq. ^5), is 
proportional to v\ cot A. Albeit the actual meridional flow profile 
inside the stars is unknown, the increase and decrease of the lati- 
tudinal flow velocity from the pole down the equator is to yield a 
peak at intermediate latitudes. In combination with the factor cot A, 
which continuously increases toward the equator, the peak of the 
density contrast will be located at low to mid latitudes. 



APPENDIX C: EVOLUTION OF SINGLE BIPOLES 

The surface flux transport simulations in Sect. l4.3l show that the in- 
termediate range of originating latitudes (Case I, cf . Table0 causes 
a higher degree of intermingling than the two cases S and L. To in- 
vestigate this aspect in more detail, we analyse the relevant mech- 
anisms in the simplified case of a single bipole evolving across 
the surface under the effect of meridional flow, differential rota- 
tion and supergranular diffusion. For a meridional flow velocity of 
vm = 100 m ■ s _1 , the butterfly diagrams in Case S, I, and L have 
the mean latitude of bipole emergence A e = 21°, 39°, and 45°, 
respectively. In each of the three cases, a single bipole is inserted 
onto a (magnetically empty) stellar surface at the mean latitude of 
emergence. The initial tilt angle is A e /2, whereas the initial total 
magnetic flux is in all cases the same (1.5 ■ 10 23 Mx). The temporal 
evolution of the polar flux density and ratio of opposite polarities 
is shown Fig. lCll The higher the (mean) latitude of emergence, the 
sooner the bipole reaches the polar region, and the stronger is the 
remaining polar flux, due to the shorter time span available for flux 
annihilation as it is transported to the poles. The trailing flux region 
(here, of positive polarity) enters the polar region before the lead- 
ing region (of negative polarity); when the latitude of emergence 
is small the timing when both polarities enter the polar region be- 
comes similar. In contrast, the higher the latitude of emergence the 
less intermingling between magnetic flux of opposite polarity is 
obtained. As the entire flux is pushed into the polar region, the flux 
ratio eventually approaches unity. Using the area under each curve 
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as a measure for the degree of intermingling, the results for the two 
lower mean latitudes of emergence are similar, whereas the degree 
of intermingling for A e = 45° is overall considerably lower. 

Since each bipole obeys Joy's law, the trailing flux region is lo- 
cated at somewhat higher latitudes than the preceding flux region of 
the opposite polarity. The meridional circulation then pushes both 
magnetic flux regions in this form toward the pole. However, dur- 
ing the poleward transport the differential rotation however rotates 
the bipole. This causes more flux of each polarity to be located at a 
common latitude, which is a requirement for intermingling to occur 
once the bipole is pushed into the polar region. If a bipole emerges 
at high latitudes, then the differential rotation has not enough time 
to act before its entry into the polar region. Hence a higher degree 
of intermingling is obtained if the initial latitude of emergence is 
low. This is in contrast to the polar flux density, which is found 
to decrease with the initial latitudes of emergence. Consequently, 
there is a particular range of latitudes of emergence and poleward 
meridional flows, which results in both high degrees of intermin- 
gling and strong polar flux densities. Optimal values occur when 
the timescales of both meridional flow and differential rotation are 
similar. 



